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Abstract:  Photodiode  nonlinearity  measurements  using  one-,  two-  and 
three -tone  measurement  systems  are  compared  with  each  other,  to 
investigate  the  comparison  accuracy  between  setups.  The  mathematical 
relationship  between  each  setup  is  analyzed,  and  data  on  multiple  devices 
are  compared  to  find  under  which  conditions  the  measurements  are 
comparable.  It  is  shown  that  the  three  measurement  systems  can  be  used 
interchangeably  only  when  the  distortion  adheres  to  the  expected 
mathematical  slopes. 
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1.  Introduction 

Photodiode  linearity  is  an  important  metric  for  analog  optical  links  requiring  high  dynamic 
range.  Nonlinearity  measurement  in  photodiodes  has  been  accomplished  using  a  number  of 
setups  usually  involving  either  heterodyning  techniques  or  Mach-Zehnder  modulators  (MZM) 
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to  create  fundamental  RF  tones  modulated  on  the  optical  input  [1-9].  Harmonic  distortion  is 
measured  using  a  pure  RF-modulated  optical  signal  incident  on  the  photodiode,  which 
generates  photodiode-induced  harmonics  in  addition  to  a  fundamental  tone.  The  measured 
powers  of  the  fundamental  and  harmonic  tones  can  then  be  extrapolated  as  a  function  of  input 
RF  power  to  calculate  the  second-  and  third-order  harmonic  output  intercept  points  (OIP2lt 
and  OIP3i,  respectively).  Intermodulation  distortion  is  measured  using  multiple  tones  to 
induce  mixing  at  the  sum  and  difference  frequencies  of  the  fundamental  tones.  Again,  OIP 
values  can  be  determined  by  extrapolating  the  measured  RF  powers  of  the  distortion  and 
fundamental. 

Previous  work  has  detailed  theoretical  mathematical  relationships  between  different  setups 
[3,5,9],  however  there  has  been  no  study  to  determine  when  these  relationships  are  physically 
true  in  photodiode  operation.  The  mathematical  relationships  make  various  assumptions  (i.e. 
memoryless,  well  behaved  nonlinearities  operating  under  sinusoidal  input)  that  may  not 
always  hold  true  in  every  case  for  every  device.  The  nonlinearities  generated  from  harmonic 
and  intermodulation  distortion  are  still  largely  not  understood-  numerous  authors  have 
presented  arguments  for  various  possible  causes  [4,7,10],  but  little  has  been  done  to  explore 
the  consistency  between  various  measurement  setups.  The  aim  of  this  work  is  to  present  three 
setups  used  to  measure  photodiode  distortion,  develop  the  mathematical  relationships  that 
allow  the  comparison  of  OIP  across  all  measurement  setups  and  measure  multiple 
photodiodes  to  determine  when  the  OIP  is  consistent  between  setups. 

With  a  multitude  of  setups  that  are  used  to  measure  nonlinearities,  there  has  been  no 
comprehensive  study  of  possible  discrepancies  between  each  measurement  system  and  the 
conditions  under  which  each  system  accurately  measures  photodiode  distortion.  Considering 
the  two-tone  MZM  setup,  Draa  et  al.  showed  that  distortion  caused  by  the  MZMs  could 
introduce  additional  nonlinearities  into  the  measurement,  making  it  inaccurate  when  using  the 
traditional  setup  [2],  Ramaswamy  et  al.  proposed  using  a  three -tone  MZM  setup  to  eliminate 
the  contribution  of  modulator-induced  distortion  to  the  measured  third  order  intermodulation 
distortion,  using  mathematical  conversions  to  relate  the  measured  data  to  the  two-tone  setup. 
The  three-tone  setup  provided  higher  OIP3  results  under  certain  conditions,  leading  the 
authors  to  conclude  that  the  three -tone  setup  is  the  measurement  system  with  less  chance  of 
MZM  nonlinearities  directly  interacting  with  the  photodiode  nonlinearities  at  high  OIP3  levels 
[9].  However,  the  accuracy  of  the  comparison  relies  on  the  mathematical  relationship  between 
the  two-  and  three -tone  measurements.  If  this  relationship  cannot  be  experimentally  verified, 
the  results  of  the  two  setups  cannot  be  compared.  Additionally,  Pan  et  al.  measured  devices 
using  both  the  two-tone  and  three-tone  MZM  setups  and  showed  little  measurement 
differences  between  them  with  their  device  even  at  high  OIP3  values  (~50dBm)  [4].  The 
heterodyning  technique  has  been  touted  for  its  lack  of  additionally  contributed  nonlinearities, 
but  is  generally  considered  to  be  the  most  complex  of  the  setups,  and  has  never  been  directly 
compared  to  modulator  measurements  [6]. 

In  all  the  works  mentioned,  there  is  no  consensus  as  to  which  measurement  system 
provides  the  most  accurate  results,  and  whether  the  data  from  one  setup  can  be  equivalently 
measured  using  another.  The  purpose  of  this  work  is  to  measure  photodiode  distortion  using 
the  three  widely  used  setups  at  a  range  of  bias  voltages  and  modulation  depths;  the  three -tone 
MZM  and  the  one-  and  two-tone  heterodyne  setups  will  be  used  to  measure  photodiode 
distortion  in  order  to  compare  nonlinearity  data  from  each  setup.  The  results  will  determine 
under  which  circumstances  the  theoretical  mathematical  relationships  can  be  applied  to  each 
setup. 

2.  Measurement  setups 

The  five -laser  heterodyne  setup  can  be  seen  in  Fig.  1.  There  are  two  pairs  of  Lightwave 
Electronics  LWE-125  diode  pumped  Nd:YAG  lasers  operating  at  a  1319nm  optical 
wavelength.  The  optical  offset  frequency  for  each  pair  is  stabilized  by  mixing  the  resulting  RF 
intermodulation  tone  (/,|/3-/,2/4)  with  an  RF  reference  signal  and  running  the  mixed  product 
through  a  phase-locked  loop  (PLL)  amplifier  circuit  back  into  the  piezo -electric  wavelength 
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control  for  one  of  the  lasers.  Each  output  is  fed  through  an  optical  attenuator  and  then 
combined  using  a  50/50  optical  coupler  and  combined  with  a  5th  laser  (Nd:YAG)  used  to 
control  the  average  photodiode  (PD)  current.  The  device  under  test  (DUT)  is  illuminated  with 
the  resulting  signal  and  the  RE  output  power  is  measured  by  an  electrical  spectrum  analyzer 
(ESA).  The  modulation  index  is  changed  by  attenuating  each  pair  of  heterodyned  lasers,  and 
adjusting  the  attenuation  of  the  5th  laser  (a5)  to  compensate  so  that  the  total  photocurrent  is 
kept  constant.  In  this  setup,  the  wavelengths  of  the  lasers  must  be  carefully  controlled  so  that 
all  unwanted  optical  intermodulation  tones  in  the  RF  frequency  domain  fall  outside  the 
bandwidth  of  the  photodiode  under  test.  The  tuning  range  of  the  LWE-125  laser  is  about 
40GHz,  so  each  set  of  heterodyned  lasers  and  the  DC  laser  were  each  tuned  to  be  20GHz  from 
each  other,  far  outside  the  3GHz  bandwidth  of  the  detector.  The  signal  produced  is  stable 
enough  that  measurements  can  be  taken  with  a  10Hz  resolution  bandwidth  (RBW),  and 
averaging  the  results  over  a  30  second  span  with  no  movement  of  the  fundamental  or 
harmonics.  The  single-tone  setup  is  configured  by  completely  tinning  off  one  pair  of 
heterodyned  lasers. 

Figure  2  shows  the  three-tone  MZM  setup.  Three  Lightwave  Electronics  model-125  diode 
pumped  Nd:YAG  lasers  are  modulated  by  three  MZMs  and  fed  through  optical  attenuators. 
The  first  two  lasers  are  combined  using  a  50/50  coupler,  and  then  to  the  third  with  another 
50/50  coupler.  Finally,  an  unmodulated  Nd:YAG  laser  is  input  to  an  optical  attenuator  and 
combined  with  the  signal  using  a  50/50  coupler.  The  output  is  fed  into  the  DUT,  and  the  RF 
output  is  measured  by  an  ESA  (Agilent  E4440A)  using  a  10Hz  RBW  and  VBW  with  a  lKHz 
span.  The  noise  floor  is  -140dBm  with  averaging,  which  gives  a  minimum  detectable 
nonlinearity  of  -120dBm  because  of  a  20dB  RF  pad  in  front  of  the  ESA.  The  MZM  setup 
utilizes  similar  laser  wavelengths  as  the  heterodyne  setup  to  remove  measurement  differences 
resulting  from  optical  wavelength  dependencies  of  the  photodiode  under  test.  For  all 
measurements,  20dB  attenuation  was  required  between  the  bias  tee  and  the  ESA  to  insure  that 
the  distortion  produced  in  the  ESA  was  below  that  of  the  PD. 


Fig.  1 .  Five  laser  two-tone  heterodyne  setup 


Fig.  2.  Four  laser  three-tone  MZM  measurement  setup 

Both  of  these  setups  were  designed  to  minimize  the  contribution  of  any  distortion 
produced  by  the  measurement  system  on  the  measured  photodiode  nonlinearity.  In  the  five 
laser  two-tone  heterodyne  setup,  additional  tones  will  be  produced  at  Xy 3  +  X2/4,  but  this  value 
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will  be  far  above  the  bandwidth  of  the  detector.  In  the  four  laser  three-tone  MZM  setup, 
second  and  third  order  distortion  terms  will  be  produced  which  have  only  FI,  F2  and  F3  terms 
contributing,  with  no  multiples  of  the  fundamental,  such  as  2-F1-F2.  This  removes  the  chance 
of  any  second  harmonics  of  the  fundamentals  contributing  to  a  third  order  distortion  product. 
In  both  setups,  we  can  be  certain  that  nonlinearity  contributed  by  the  measurement  system  is 
minimized. 

3.  Mathematical  relationships 

To  determine  the  mathematical  relationship  between  the  data  from  each  setup,  a  generalized 
photodiode  responsivity  nonlinearity,  55,  will  be  designated.  We  can  express  photocurrent  as  a 
function  of  optical  power  by: 


i  =  pumx) 


(i) 


We  will  define  9I(x)  as  a  nonlinear  transfer  function  of  a  physical  variable  x  that  can  be 
approximated  using  a  Taylor  series  around  some  arbitrary  value,  a,  up  to  the  second  order: 


9f(x)  «  T?(a)  +  9f'(fl)(x-fl)+  ^  ^a\x-a)2 


(2) 


The  incident  optical  power  in  each  setup  will  be  defined  by  a  sum  of  sinusoids. 


P(n  =  p 


(  k 

M 

i  + 

’’Tin,  sin (®;?) 

V  /=! 

J3 

(3) 


where  k  =  1  for  the  one-tone  heterodyne,  k  =  2  for  the  two-tone  heterodyne  and  k  =  3  for  the 
three -tone  MZM  setup,  P  is  the  average  power,  ?n,  is  the  modulation  index,  and  ®(  is  the 

RF  tone  frequency.  In  the  heterodyne  setup,  ml  can  be  determined  by  simply  measuring  the 
average  power  contributions  from  each  of  the  heterodyned  laser  pairs  and  the  DC  laser.  In  the 
modulator  setup,  m,  can  be  calculated  using  the  RF  input  power  to  the  modulator,  optical 
input  power,  and  modulator  transfer  function  [11],  We  can  assume  x  is  some  slowly  varying 
function  of  P(t)  and  that  small  signal  conditions  apply  so  that: 


x«d-P(f)  (4) 

where  S  is  some  constant.  Substituting  Eq.  (4)  and  Eq.  (2)  into  Eq.  (1)  we  have: 

/  =  ctP(r)  +  c2P2  (?) + c3P 3  (?)  (5) 


where  c,  =  9f(a)  -9T(a)a  + 


SR'(a)  2  ,9T(«)  *  ,  SR'(a)  C2  ^ 

-a  ,  c2  =  :H  (a)d  -2 - da,  and  c3  =  — r — d  .  For  a 


2  2  2 
one -tone  signal  where  1  =  1,  we  can  substitute  the  expanded  form  of  Eq.  (3)  into  Eq.  (5). 


I  =  paAc\  [1  +  '«i  sinfOj? )]  +  c2 


1  +  2  m,  sin(®1?)  +  ^-(l-cos^®,?)) 


3  iy[~ 

1  +  3mt  sinfOj?)  +  — (1  - 008(2®,?))  +  (3 sin(®,?)  - 8111(3®,?)) 


(6) 


From  Eq.  (6),  the  coefficients  for  the  fundamental,  second  order  and  third  order  components 
can  be  extracted.  Repeating  this  series  of  steps  for  k  =  2  and  3 ,  the  resulting  intermodulation 
and  harmonic  terms  in  each  measurement  system  can  be  isolated  and  compared  with  each 
other.  This  small  signal  approximation  comparison  holds  for  any  well-behaved  single 
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photodiode  nonlinearity,  so  that  for  every  dB  increase  in  the  fundamental,  the  nth  order 
distortion  increases  by  n  dB. 

When  the  one-,  two-  and  three-tone  cases  are  compared  with  each  other,  similar  ordered 
distortion  in  each  case  can  be  compared  while  keeping  the  modulation  index  constant.  The 
second  order  relationships  using  the  two-tone  setup  as  the  basis,  are  described  by 

HD2U  +  6  =  IMD22i  =  IMD23l  (dBm)  (7) 

where  HD2U  is  the  second  order  harmonic  distortion  for  the  one-tone  setup  and  IMD22t  and 
IMD23i  are  the  second  order  intermodulation  distortion  for  the  two-tone  and  three-tone  setups 
respectively  [12].  For  the  third  order  distortion,  the  relationships  using  the  two -tone  setup  as 
the  base,  are  described  by 

HD3h  +  9.54  =  mm2t  =  IMD3  .,  -  6  (dBm)  (8) 

where  HD3U ,  IMD3 and  IMD3it  are  the  third  order  distortion  for  the  one-,  two-  and  three- 
tone  setups  respectively  [12].  Finally,  the  relationship  between  OIP2  and  OIP3  for  each  setup 
can  be  determined  by. 


OIP2  =  2 F-  IMD2  (dBm)  (9a) 

OIP3  =  |f-^ IMD3  (dBm)  (9b) 


where  F  is  the  average  power  of  the  fundamentals,  in  dBm,  and  IMD2I IMD3  can  also  be 
the  harmonic  distortion  ( HD2I HD3)  of  the  one-tone  setup.  Using  Eq.  (9a),  the  relationship 
between  OIP2,  using  the  two-tone  setup  as  the  base  is 


OIP2u  -  6  =  OIP2lt  =  OIP23i  (dBm)  ( 1 0) 

Similarly,  using  Eq.  (9b)  the  relationship  between  OIP3  using  the  two-tone  setup  as  the 
base,  in  dBm  units  of  power,  is 

OIP3h  -A. 11  =  OIP32i  =  OIP3.,  +  3  (dBm)  (11) 

Using  Eq.  (10)  and  Eq.  (11),  the  data  from  the  one-  and  three-tone  measurement  systems  can 
be  shifted  according  to  the  constants  so  that  comparisons  between  the  three  systems  can  be 
made. 

4.  Results  and  analysis 

4.1  Applied  Electronics  PD3000 

A  3GHz-bandwidth  photodiode  (Applied  Electronics  model:  PD3000)  was  measured.  The 
measurement  was  taken  around  a  100MHz  center  frequency  with  10MHz  spacing  for  the  two- 
and  three -tone  setups,  at  10mA  average  photocurrent.  The  tones  were  set  to  slightly  offset 
frequencies  (i.e.  89.59MHz  instead  of  90.00MHz)  instead  of  integers  to  avoid  any  extraneous 
signals  present  in  the  environment  that  might  mix  and  add  to  distortion  at  even  numbers  (i.e. 
10MHz,  100MHz).  All  second  and  third  order  distortion  products  were  measured  as  a  function 
of  photodiode  bias  voltage  and  modulation  depth.  In  general  practice,  OIP  is  extracted  by 
plotting  the  fundamental  and  second  and  third  order  distortion  as  a  function  of  modulation 
depth  or  input  power  and  then  fitting  each  with  a  slope  of  one,  two,  and  three  respectively 
according  the  mathematical  relationships  established  in  section  3.  Assuming  this  is  the  case, 
local  OIP2  and  OIP3  were  calculated  at  each  bias  voltage  point  using  Eq.  (9a)  and  Eq.  (9b)  for 
the  same  modulation  depth  across  each  setup.  The  data  were  plotted  at  the  same  input 
modulation  depth  across  all  setups  to  get  the  most  accurate  local  OIP2  and  OIP3  comparison 
between  setups. 
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Figure  3  shows  the  OIP2  versus  bias  voltage  for  all  second  order  distortion  in  the  three 
setups  separated  into  areas  A  and  B.  The  plot  contains  the  second  order  intermodulation 
distortions  that  add  (i.e.  FI  +  F2)  as  well  as  those  that  subtract  (i.e.  F2-F1).  The  one-tone 
heterodyne  setup  provides  the  second  order  harmonic  (2- FI)  at  200MHz.  To  compare  the 
three  measurement  setups,  the  01P2h  and  OIP 2it  data  are  shifted  according  to  Eq.  (10)  to 
()IP2ni .  The  three  measurement  setups  are  designated  by  color,  red  for  the  one-tone,  blue  for 
the  two-tone,  and  green  for  the  three-tone.  The  color  scheme  for  each  setup  will  be  maintained 
throughout  the  remaining  figures.  At  low  bias  voltages  (area  A),  the  data  clearly  line  up  for 
low  frequency  (solid  points)  and  high  frequency  (open  points)  distortion.  At  high  bias 
voltages  (area  B),  the  data  at  high  frequencies  for  the  three-tone  setup  remain  consistent. 
There  is  a  large  divergence  between  all  three  setups  from  5  to  6V.  For  the  low  frequency 
distortion  at  high  voltages,  F3-F1  and  F3-F2  for  the  three-tone  setup  merge  together  along 
with  F2-F1  of  the  two-tone,  where  before  F3-F1  was  significantly  lower  than  the  distortion 
closer  to  10MHz.  The  difference  between  the  low  and  high  frequencies  will  be  addressed 
later. 


Bias  Voltage  (V) 

Fig.  3.  OIP2  versus  bias  voltage  for  a  PD3000  photodiode,  using  one-,  two-  and  three-tone 
measurement  setups.  Fundamental  frequency  f  =  lOOMFlz,  with  frequency  spacing  Af  = 
lOMFlz  and  total  photocurrent  =  10mA.  One-  and  three-tone  data  are  shifted  according  to  Eq. 

(10)  to  mathematically  match  the  two-tone  data.  One-tone  data  for  5MHz  fundamental  included 
(red  solid  triangles). 

In  Fig.  4,  OIP3  is  plotted  versus  bias  voltage  from  0  to  10V  for  intermods  with  frequencies 
close  to  the  fundamental  frequency.  In  this  instance  the  MZM  setup  contributes  three  sets  of 
data,  FI  +  F2-F3,  F1-F2  +  F3,  and  F2  +  F3-F1,  while  the  two-tone  heterodyne  setup 
contributes  data  for  2-F1-F2  and  2-F2-F1.  To  compare  the  data  sets  from  each  measurement 
setup,  the  one-  and  three-tone  data  are  shifted  as  per  Eq.  (11),  to  compare  them  to  the  two- 
tone  setup.  In  Fig.  4,  the  data  match  up  well  in  area  A  (  +  ldB).  In  area  B,  the  two-tone  setup 
and  three-tone  setup  no  longer  overlap.  Also,  the  different  OIP3  values  within  each  setup 
begin  to  deviate  from  each  other.  For  the  three-tone  setup  there  is  about  ldB  between  each 
data  set,  with  the  OIP3  decreasing  as  the  frequency  of  the  distortion  increases  (78MHz, 
102MHz,  and  122MHz).  For  the  two-tone  setup,  the  data  cross  back  and  forth  with  no  clear 
pattern.  The  largest  discrepancy  between  the  measurement  systems  is  8dB  corresponding  to  a 
16dB  difference  in  distortion,  which  is  quite  substantial. 

Similarly,  the  OIP3  data  for  intermods  close  to  the  third  harmonic  frequency  are  plotted  in 
Fig.  5.  The  three-tone  setup  provided  FI  +  F2  +  F3,  the  two-tone  FI  +  2-F2  and  2- FI  +  F2, 
and  the  one-tone  3F1.  The  data  match  well  in  area  A  as  in  Fig.  4;  however  the  three  setups 
begin  to  diverge  at  higher  bias  voltages  in  area  B  as  seen  before  in  Figs.  3  and  4.  After  5V 
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there  is  a  large  divergence  from  any  pattern  in  the  data.  The  one -tone  data  in  particular  are 
much  lower  than  the  other  two  setups. 


Fig.  4.  OIP3  versus  bias  voltage  for  a  PD3000  photodiode,  using  one-,  two-  and  three-tone 
measurement  setups.  Fundamental  frequency  f  =  100MHz,  with  frequency  spacing  Af  = 
10MHz  and  total  photocurrent  =  10mA.  Only  distortion  near  f  is  plotted.  One-  and  three-tone 
data  are  according  to  Eq.  (11)  to  mathematically  match  the  two-tone  data.  One-tone  data  for 
33MHz  fundamental  included  (red  solid  triangles). 


Fig.  5.  OIP3  versus  bias  voltage  for  a  PD3000  photodiode,  using  one-,  two-  and  three-tone 
measurement  setups.  Fundamental  frequency  f  =  100MHz,  with  frequency  spacing  Af  = 

10MHz  and  total  photocurrent  =  10mA.  Only  distortion  near  3-f  is  plotted.  One-  and  three-tone 
data  are  according  to  Eq.  (1 1)  to  mathematically  match  the  two-tone  data. 

The  OIP3  data  are  plotted  separately  for  low  (near  the  fundamental)  and  high  frequency 
cases  (near  3  times  the  fundamental)  because  the  shapes  of  the  curves  do  not  match  up  at  the 
two  sets  of  frequencies  (similar  to  the  OIP2  case).  In  the  second  order  case,  the  low  frequency 
OIP2s  (around  10MHz)  are  significantly  higher  than  the  high  frequency  OIP2s.  The  OIP2  that 
results  in  double  the  frequency  spacing  (~20MHz)  is  slightly  lower  than  the  10MHz  OIP2s.  In 
general,  the  low  frequency  cases  have  less  distortion  (and  therefore  a  higher  OIP)  than  the 
high  frequency  cases.  This  frequency  dependence  appears  to  be  more  pronounced  at  low 
frequencies  (around  10MHz)  than  at  high  frequencies  (around  200MHz).  Therefore,  there  is  a 
significant  difference  in  OIP2  between  F2-Fl(10MHz)  and  F3-F1  (20MHz)  for  the  three-tone 
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case,  with  a  separation  of  10MHz.  With  similar  separation  at  high  frequencies  there  is  no 
difference  in  OIP2  for  FI  +  F2  (191MHz)  and  FI  +  F3  (203MHz).  The  disparity  between 
same  order  distortion,  as  a  function  of  frequency  leads  to  the  conclusion  there  is  some 
frequency  dependence  within  a  certain  order  of  photodiode  nonlinearity. 

Additional  data  are  plotted  in  Figs.  3  and  4  to  analyze  this  frequency  dependence.  The 
OIP2  from  the  one-tone  setup  (red  solid  triangles)  with  a  5MHz  fundamental  resulting  in  a 
second  order  harmonic  at  10MHz  is  plotted  in  Fig.  3.  In  area  A,  the  data  line  up  well  with  the 
low  frequency  data  from  the  two-  and  three-tone  setups,  despite  originating  from  a  much 
lower  fundamental  frequency.  Data  are  plotted  in  Fig.  4  from  the  one -tone  setup  (red  solid 
triangles)  with  a  fundamental  of  33MHz,  which  results  in  a  third  order  harmonic  at  99MHz. 
The  data  match  well  with  the  other  nonlinearities  around  100MHz  in  area  A.  Again,  in  both 
cases,  the  harmonic  distortion  deviates  at  higher  bias  voltages  in  area  B.  The  general 
patternless  deviation  of  all  distortion  at  higher  bias  voltage  is  probably  a  result  of  nonlinearity 
characteristics  specific  to  this  photodiode,  since  we  also  observed  photodiodes  where  the 
distortion  matches  up  well  at  high  bias  voltages. 


Fig.  6.  Plots  of  the  electric  field  vs.  time  for  one-tone  (top),  two-tone  (middle)  and  three-tone 
(bottom)  measurements. 


Fig.  7.  Fundamental  and  selected  second  and  third  order  distortion  power  versus  modulation 
depth  at  2.5V  bias  for  a  PD3000  photodiode.  Fundamental  frequency  f  =  100MHz,  with 
frequency  spacing  Af  =  10MHz  and  total  photocurrent  =  10mA.  One-  and  three-tone  data  are 
according  to  Eq.  (10)  and  Eq.  (1 1)  to  mathematically  match  the  two-tone  data. 

The  large  difference  in  low  and  high  frequency  second  order  distortion  values  can  be 
attributed  to  the  differences  in  electric  field  that  the  diode  sees  between  the  measurement 
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systems.  The  plots  of  the  electric  field  for  the  one-tone,  two-tone  and  three-tone  sinusoidal 
measurement  systems  are  shown  in  Fig.  6.  The  data  in  Fig.  3  and  Fig.  4  from  the  one-tone 
setup  at  different  fundamentals  (5MHz  and  33MHz)  demonstrate  that  the  OIP2  and  OIP3  are 
not  necessarily  a  strong  function  of  the  fundamental  frequencies.  Instead  the  nonlinearity  may 
be  a  function  of  the  time  dependent  envelopes  of  the  electric  field  rather  than  from  the  higher 
frequency.  Since  low  frequency  nonlinearity  is  dependent  on  the  envelope  frequency,  and  the 
high  frequency  nonlinearity  is  dependent  on  the  internal  “carrier”  frequency,  the  resulting 
distortion  levels  will  differ  to  reflect  the  large  difference  in  frequency  they  originate  from. 


(a)  (b) 

Fig.  8.  Slope  of  (a)  IMD2/HD2  and  (b)  IMD3/HD3  vs.  bias  voltage  for  PD3000  for  all  three 

setups. 

Now  that  the  diode  has  been  analyzed  over  a  range  of  bias  voltages  to  determine  the 
correlation  of  the  setups,  the  relationship  of  the  output  RF  power  and  modulation  depth  can  be 
discussed.  In  Fig.  7,  fundamental,  second  and  third  order  distortion  power  is  plotted  versus 
modulation  depth  at  2.5V  bias,  corresponding  to  area  A  in  Figs.  3-5,  where  the  data  from  the 
different  setups  were  well-correlated.  Again,  the  one-  and  three-tone  distortion  second  and 
third  order  distortion  data  have  been  shifted  to  the  two-tone  setup,  as  per  Eq.  (7)  and  Eq.  (8) 
respectively.  In  the  plot,  the  second  order  data  clearly  line  up  for  low  and  high  frequencies 
over  the  entire  measured  range  of  modulation  depths.  Additionally  the  data  are  all  within  3% 
deviation  of  the  expected  slope  value  of  two.  The  slopes  of  the  3rd  order  distortion  are 
similarly  behaved  (slope  of  three)  across  all  measurement  setups.  All  of  the  data  agree  to 
within  +  ldB  at  different  modulation  depths,  indicating  that  all  setups  are  measuring  the  same 
nonlinearity,  and  therefore  producing  the  same  local  OIP.  In  this  case,  we  conclude  that  the 
measurement  setups  can  be  used  interchangeably  since  the  mathematical  relationships  in  Eq. 
(10)  and  Eq.  (11)  hold  true.  To  illustrate  the  desired  measurement  behavior  more  clearly,  the 
slopes  of  the  measured  IMD2/HD2  and  IMD3/HD3  are  plotted  vs.  bias  voltage  in  Fig.  8a  and 
Fig.  8b  respectively.  The  slopes  match  up  to  2  very  well  up  in  area  A  in  Fig.  8a.  In  Fig.  8b 
there  is  some  discrepancy  at  0.5V,  IV  and  3.5V,  but  the  general  trend  does  not  diverge  until 
higher  voltages  in  area  B.  The  larger  error  in  the  slope  values  for  the  third  order  distortion  is 
attributed  to  fact  that  the  third  order  distortion  is  a  more  complex  process  than  the  second 
order,  and  thus  is  not  as  consistent  with  the  ideal  theoretical  behavior.  When  the  slopes 
diverge  significantly  from  their  expected  values,  OIP  becomes  impossible  to  calculate,  and 
Eq.  (10)  and  Eq.  (11)  are  no  longer  useful  in  comparing  the  three  setups.  For  Figs.  3-5  this 
makes  the  OIP  values  for  bias  voltage  higher  than  4V  invalid. 
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4.2  Additional  Devices 


In  order  to  confirm  that  the  behavior  in  4.1  is  not  a  singular  case  for  the  diode  tested,  two 
other  devices  were  measured  using  the  three  setups.  A  second  commercial  diode  from  Applied 
Electronics  (Model:  PD  1000)  was  tested  at  10mA  at  a  center  frequency  of  100MHz  with 
10MHz  spacing  from  0  to  10V  bias.  The  PD1000  is  similar  to  the  PD3000,  but  has  a 
bandwidth  of  1GHz  instead  of  3GHz.  The  OIP2  vs.  bias  voltage  is  plotted  in  Fig.  9a.  The 
OIP2  is  similar  to  that  of  the  PD3000  in  Figs.  4  and  5,  where  the  setups  match  up  well  in  area 
A.  The  slope  of  the  distortion  is  plotted  in  Fig.  9b  where  the  slope  is  close  to  2  in  area  A,  but 
exhibits  more  discrepancy  than  in  the  PD3000  case,  particularly  at  0.5  V,  where  all  frequencies 
show  a  spike  in  slope.  The  OIP3  for  the  PD1000  was  also  measured  and  exhibited  similar 
behavior  to  the  PD3000  and  is  thus  not  graphed.  As  before  in  area  B  the  OIP2  diverges, 
meaning  the  OIP2  cannot  be  extracted  as  the  slopes  do  not  match  the  mathematical  values. 


Fig.  9.  (a)  OIP2  vs.  bias  voltage  for  PD1000  at  center  frequency  100MHz  with  10MHz  spacing 
at  10mA  photocurrent  for  one-,  two-  and  three-tone  setups  and  (b)  IMD2/H23  slope  vs.  bias 
voltage  for  PD1000.  In  (a),  one-  and  three-tone  data  are  according  to  Eq.  (10)  to 
mathematically  match  the  two-tone  data. 

A  third  commercial  device  was  measured  that  was  not  from  Applied  Electronics,  which 
will  be  referred  to  as  Device  C  and  has  a  bandwidth  of  20GHz.  Device  C  was  measured  at  the 
same  conditions  as  the  PD3000  in  section  4.1.  The  OIP2  vs.  bias  voltage  is  plotted  in  Fig.  10a. 
In  this  case,  OIP2  at  all  distortion  frequencies  matches  well  in  area  B,  at  higher  bias  voltages, 
but  is  inconsistent  in  area  A,  at  lower  voltages.  To  confirm  if  the  same  conditions  used 
previously  apply,  the  slope  of  the  second  order  distortion  is  plotted  in  Fig.  10b.  It  is  clear  that 
the  slopes  are  all  approximately  2  in  area  B,  where  the  OIP2  matches  well  between  all  setups, 
where  as  in  area  A  the  slopes  do  not  match  up  and  therefore  the  OIP2  extracted  is  invalid 
according  to  the  mathematical  definition.  Since  for  Device  C  the  setups  produce  consistent 
results  at  a  different  range  of  voltages  than  for  the  PD3000  and  PD  1000,  there  is  unlikely  to 
be  a  general  explanation  for  the  discrepancy  of  distortion  slope  from  its  expected  value.  We 
can  conclude  that  nonlinearity  measurements  made  on  different  setups  cannot  be  cross 
compared  unless  the  distortion  produces  the  correct  theoretical  slope  vs.  modulation  depth, 
and  in  order  to  verify  the  slope,  detailed  distortion  measurements  as  a  function  of  modulation 
depth  must  be  made.  Deviation  of  the  distortion  slope  from  the  expected  value  could  be 
caused  by  several  factors  in  the  photodiode.  Nonlinearities  identified  previously  in  [4,7,10], 
may  provide  higher  order  terms  that  dominate  at  a  lower  order  distortion,  and  phasing/de¬ 
phasing  between  two  competing  nonlinearities  can  cause  cancellation  [13]  in  the  distortion 
measurement.  In  the  mathematical  model,  the  assumption  was  made  that  high  order  terms 
would  not  dominate;  however  slopes  higher  than  the  ideal  value  suggest  the  higher  order 
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nonlinearities  should  be  investigated  both  mathematically,  by  expanding  the  Taylor  series,  and 
experimentally  in  the  future. 


Fig.  10.  (a)  OIP2  vs.  bias  voltage  for  Device  C  at  center  frequency  100MHz  with  10MHz 
spacing  at  10mA  photocurrent  for  one-,  two-  and  three-tone  setups  and  (b)  IMD2/HD2  slope 
vs.  bias  voltage  for  Device  C.  In  (a),  one-  and  three-tone  data  are  according  to  Eq.  (10)  to 
mathematically  match  the  two-tone  data. 


5.  Conclusion 

Nonlinearity  was  measured  (OIP2  and  OIP3)  for  the  PD3000  photodiode  using  three  setups,  a 
one -tone  heterodyne,  two-tone  heterodyne  and  three-tone  MZM  design.  The  data  from  each 
setup  were  measured  over  a  range  of  10V  bias  voltage  with  a  100MHz  center  frequency  and 
10MHz  separation  between  tones.  The  results  matched  well  between  all  setups  from  low  bias 
voltages,  where  the  diode  exhibited  the  expected  behavior  versus  modulation  depth,  verifying 
the  mathematical  relationships  developed  in  section  3.  At  higher  bias  voltages,  the  results 
diverged  for  each  setup  and  were  no  longer  correlated  to  each  other.  The  slopes  versus 
modulation  depth  no  longer  adhered  to  the  mathematical  relationships  of  two  for  OIP2  and 
three  for  OIP3,  invalidating  the  three  setups  as  consistent  measurements.  Data  were  also  taken 
using  the  one -tone  setup  demonstrating  a  nonlinearity  correlation  resulting  from  the  frequency 
at  which  the  distortion  occurs,  as  opposed  to  the  frequency  of  the  fundamental.  Additionally, 
data  were  taken  on  two  other  devices,  in  which  the  theoretical  mathematical  relationship 
requirement  confirmed  the  consistency  of  the  different  setup  measurements.  We  have  shown 
that  while  all  three  measurements  systems  are  capable  of  producing  accurate  nonlinearity 
measurements,  they  do  not  do  so  under  all  circumstances.  In  order  to  insure  the  accuracy  of 
the  measurement,  more  than  a  “spot”  check  of  OIP  at  a  single  modulation  depth  is  needed. 
Even  in  the  absence  of  setup  nonlinearity  (MZM  nonlinearity),  device  nonlinearities  such  as 
higher  order  distortion  products  or  possible  phase  cancellation  between  competing 
mechanisms  can  affect  the  measured  slope  as  a  function  of  modulation  depth,  making  the 
extraction  of  OIP  impossible  per  the  mathematical  definition.  From  our  data,  we  can  assume 
that  as  long  as  the  slopes  of  the  nonlinearities  match  their  expected  values  over  a  range  of 
modulation  depths,  the  three  setups  will  measure  the  same  nonlinearity,  and  all  of  them  can  be 
interchanged. 
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